clear
u=-3:0.01:10;
for k=1:1301
   if u(k)>0
    i(k)=1;
   else
    i(k)=0;
   end
end;
plot(u,i);
legend('Uvd / Ur');